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Abstract 

Chiral symmetry breaking may exhibit significantly different patterns in two chiral limits: Nj = 
2 massless flavours {m u = = 0, m s physical) and Nf = 3 massless flavours {m u = = m s = 0). 
Such a difference may arise due to vacuum fluctuations of ss pairs related to the violation of the 
Zweig rule in the scalar sector, and could yield a numerical competition between contributions 
counted as leading and next-to-leading order in the chiral expansions of observables. We recall 
and extend Resummed Chiral Perturbation Theory (RexPT), a framework that we introduced 
previously to deal with such instabilities: it requires a more careful definition of the relevant ob- 
servables and their one-loop chiral expansions. We analyse the amplitudes for low-energy tttt and 
nK scatterings within RexPT, which we match in subthreshold regions with dispersive representa- 
tions obtained from the solutions of Roy and Roy-Steiner equations. Using a frequentist approach, 
we constrain the quark mass ratio as well as the quark condensate and the pseudoscalar decay 
constant in the Nj = 3 chiral limit. The results mildly favour significant contributions of vacuum 
fluctuations suppressing the Nf = 3 quark condensate compared to its Nf = 2 counterpart. 



1 Work supported in part by the EU Contract No. MRTN-CT-2006-035482, "FLAVIAnet" . 
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1 Introduction 



A striking feature of the Standard Model consists in the mass hierarchy obeyed by the light quarks: 

m u ~ m d <C m s ~ A QC d < A H , (1) 

where Aqqd is the characteristic scale describing the running of the QCD effective coupling and 
Ah ~ 1 GeV the mass scale of the bound states not protected by chiral symmetry. Therefore, the 
strange quark may play a special role in the low-energy dynamics of QCD: 

i) it is light enough to allow for a combined expansion of observables in powers of m u ,md,m s 
around the Nf = 3 chiral limit (meaning 3 massless flavours): 

Nf = 3 : m u = m d = m s = , (2) 

ii) it is sufficiently heavy to induce significant changes in order parameters from the Nf = 3 
chiral limit to the Nf = 2 chiral limit (meaning 2 massless flavours): 

Nf = 2 : m u = = m s physical , (3) 



Hi) it is too light to suppress efficiently loop effects of massive ss pairs (contrary to c, b, t quarks). 

These three arguments suggest that ss sea pairs may play a significant role in chiral dynamics, 
leading to different patterns of chiral symmetry breaking in Nf = 2 and Nf = 3 chiral limits. Then, 
chiral order parameters such as the quark condensate and the pseudoscalar decay constant: 

S(iV / ) = -lim(^), F 2 (Nf)=l\mF*, (4) 

Nf Nf 

would have significantly different values in the two chiral limits (lim^Vj denoting the chiral limit 
with Nf massless flavours). 

The role of ss-pairs in the structure of QCD vacuum is a typical loop effect : it should be 
suppressed in the large- iV c limit, and it can be significant only if the Zweig rule is badly violated in 
the vacuum (scalar) channel J PC = ++ . On general theoretical grounds pQ, one expects ss sea- 
quark pairs to have a paramagnetic effect on chiral order parameters, so that they should decrease 
when the strange quark mass is sent to zero : for instance, S(2;m s ) > E(2;m s = 0), and similarly 
for F 2 , which translates into the paramagnetic inequalities: 

£(2)>E(3), F 2 (2)>F 2 (3). (5) 

However, the size of this paramagnetic suppression is not predicted. Thus, it is highly desirable to 
extract the size of the chiral order parameters in Nf = 2 and Nf = 3 limits from experiment. 

Recent data on tttt scattering [2j together with older data and numerical solutions of the Roy 
equations [3] allowed us to determine the Nf = 2 order parameters expressed in suitable physical 
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units 



X(2) 



(m u + m d )Y,(2) 



= 0.81 ± 0.07 



(6) 



Z(2) 




TT 



= 0.89 ±0.03. 



(7) 



A different analysis of the data in ref. [2], with the additional input of dispersive estimates for the 
(non-strange) scalar radius of the pion, led to a larger value of X{2) [5]. X(2) and Z{2) seem fairly 
close to 1, so that corrections related to m u ,md ^ (while m s remains at its physical value) have 
a limited impact on the low-energy behaviour of QCD. In turn, two-flavour Chiral Perturbation 
Theory (%PT) [6], which consists in an expansion in powers of m u and rrid around the Nf = 2 
chiral limit, would not suffer from severe problems of convergence q 

Two-flavour xPT [S] deals only with dynamical pions in a very limited range of energy. In 
order to include K- and ^-mesons dynamically and extend the energy range of interest, one must 
use three-flavour \PT [9] where the expansion in the three light-quark masses starts around the 
Nf = 3 vacuum m u = = m s = 0. From the above discussion, large vacuum fluctuations of 
ss pairs would have a dramatic effect on Nf = 3 chiral expansions. The leading-order (LO) term, 
which depends on the 0(p 2 ) low-energy constants F 2 (S) and S(3), would be damped. On the 
other hand, next-to-leading-order (NLO) corrections could be enhanced, in particular those related 
to Zweig-rule violation in the scalar sector. For instance, the Gell-Mann-Oakes-Renner relation 
would not be saturated by its LO term and would receive sizable numerical contributions from 
terms counted as NLO in the chiral counting. 

Unfortunately, the experimental data on K- and 77-decays are not accurate enough to assess the 
role of ss pairs in the iVy-dependence of chiral symmetry breaking in a very precise way. However our 
understanding of ttK scattering at low energies has been improved recently through the re-analysis 
of dispersive Roy-Steiner equations [10] . A rapid analysis of its results in the framework of three- 
flavour xPT hinted at significant vacuum fluctuations encoded in some 0(p ) chiral couplings, 
which calls for a more detailed analysis of the ttK system. Interesting information can also be 
obtained from our current knowledge of ttit scattering, which we will include in our study. 

To perform such an analysis, we develop and modify the framework presented in refs. [Tj llll[T2] . 
Specifically, our work differs from ref. [12] on three points: we consider not only irir- but also ttK- 
scattering, our observables are the values of the amplitudes in unphysical regions rather than 
subtraction constants of dispersion relations, the matching between theoretical and experimental 
representations is performed in a frequentist approach, not in a Bayesian framework. 

In sec. 2, we motivate and explain Resummed Chiral Perturbation Theory (Re^PT), a frame- 
work designed to derive three-flavour chiral series at one loop, in which vacuum fluctuations of ss 
pairs are resummed. In sec. 3, we apply RexPT to tttt- and vri^-scattering amplitudes. In sec. 4, 
we explain how we determine the same amplitudes dispersively in subthreshold regions, building 
upon the solutions of Roy and Roy-Steiner equations O [10] - In sec. 5, we discuss the matching of 

2 Let us stress that new data of high accuracy are expected from the NA48/2 collaboration soon, which could 
affect these results significantly [7]. Recent lattice simulations with two-flavour dynamical quarks [8] may help to 
understand some aspects of these questions, even though the results are preliminary and rather delicate to interpret. 
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the chiral and dispersive results within a frequentist approach |13j . and in sec. 6, we present our 
results for the order parameters of Nf = 3 chiral symmetry breaking. In sec. 7, we summarise and 
discuss our results. Appendices are devoted to the expression of scattering amplitudes in Re^PT, 
their evaluation from Roy and Roy-Steiner equations and the treatment of correlated data. 

2 Resummed Chiral Perturbation Theory 

We start by describing in more detail the framework introduced in refs. [TJ \TT\ [T2] to expand 
observables around the Nf = 3 chiral limit in the case of significant vacuum fluctuations. We take 
this opportunity to extend this framework to deal with energy-dependent quantities. 

2.1 Convergence of observables 

In the introduction, we have emphasised the possibility for three-flavour chiral series to exhibit 
a rather unusual behaviour, with a numerical competition between leading and next-to-leading 
order. In ref. [12], we called instability of the expansion such a numerical competition between 
terms of different chiral counting. A naive argument based on resonance saturation suggests that 
higher orders in the chiral expansion should be suppressed by powers of (F^/Ah) 2 . However, such 
an argument does not apply to a leading-order contribution proportional to S(3) [TS]: there is no 
resonance that could saturate the quark condensate. Therefore we expect to encounter three-flavour 
chiral expansions with a good overall convergence: 

A = A LO + A NLO + A5A , 5A4Z.1, (8) 

but the numerical balance between the leading order Alo an d the next-to-leading order A/vlo 
depends on the importance of vacuum fluctuations of ss pairs. 

At the level of 0(p 4 ) Nf = 3 chiral perturbation theory, the size of the vacuum fluctuations is 
encoded in the low-energy constants (LECs) L4 and Lq whose values remain largely unknown. For 
a long time [9], they have been set to zero at an arbitrary hadronic scale (typically the r/-mass) 
assuming that the Zweig rule held in the scalar sector. More recent but indirect analyses based 
on dispersive methods [I0l El HE] suggest values of L4 and Lq which look quite modest but are 
sufficient to drive the Nf = 3 order parameters S(3) and i ?2 (3) down to half of their Nf = 2 
counterparts S(2) and -F 2 (2), leading to Alo — Aj^lo as recalled in sec. 12.31 In addition, two-loop 
analyses [HI [201 1211 H2] led to values of L4 and Lq off large- N c expectations. 

Unstable Nf = 3 chiral expansions (Alo ~ Anlo) demand a more careful treatment than in 
two-flavour xPT where such instabilities are seemingly absent. For instance, it would be wrong to 
believe that the chiral expansion of 1/A converges nicely El- This might induce the observed problems 
of convergence in recent two-loop computations [HH HSU EI] : the latter treat the fluctuations 
encoded in L4 and Lq as small and are not designed to cope with a large violation of the Zweig 
rule in the scalar sector, leading to instabilities of the chiral series. 

3 This would be equivalent to claiming that 1/(1 + a;) ~ 1 — x is a reasonable approximation for x = O(l). 
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Observables with a good convergence in the sense of eq. ([8]) form a linear space, which we 
identify with connected QCD correlators of axial/vector currents and their derivatives, away from 
kinematic singularities. This choice promotes some "good" observables that can be extracted from 
such correlators, such as Fp and FpMp (P = ir, K, rj) : LO and NLO may compete, but there 
should be only a tiny contribution from NNLO and higher orders. On the contrary, the chiral 
expansion of Mp (ratio of the former quantities) may exhibit a bad convergence. Similarly, the 
good observable associated with a form factor Fp^Q describing a transition from a pseudoscalar 
meson P to a meson Q will be FpFqFp^q, where the decay constants Fp and Fn stem from 
wave-function renormalisation factors in the LSZ reduction formula. 



2.2 One-loop bare expansion of QCD Green functions 



In a previous work [12] . we proposed a framework to deal with chiral expansions in the case of 
large fluctuations, by resumming the terms containing the Zweig-rule violating LECs L4 and Lq. 
This framework, which we will call Resummed Chiral Perturbation Theory (RexPT), includes 
consistently the alternatives of large and small vacuum fluctuations. In this section, we explain 
how to expand a good observable at one loop in RexPT. We addressed only energy-independent 
quantities in ref . [112] , where we explained in detail the similarities and differences of our approach 
with respect to Generalized Chiral Perturbation Theory [151 116] . 

We start from the one-loop generating functional for three-flavour xPT [H]: 

Z = Z t + Z u + Z A + ... (9) 

where the ellipsis stands for NNLO contributions. The three terms of the one-loop generating 
functional are: 



Zt is the sum of 0(p 2 ) and 0(p A ) tree graphs, and of tadpole contributions: 

Zt = yJd X fi-TL^o g Mp\a^ P (10) 




dx £4. i 

— J u I uvi- r 6 fj,- 1 j 

where Fq = F(3), 

Fo = F(3), B = ^ ,r = ^ (11) 
<7 A and a x collect source terms for vector / axial currents and scalar / pseudoscalar densities, 

o 2 

and £4 is the 0(p ) chiral Lagrangian with renormalised couplings L\ and HJ. Mp denotes 
the 0(p 2 ) contribution to the (squared) mass of the Goldstone boson P: 

mI= Y{2,)Ml , M K = r -^Y(3)M 2 , m1= ^±±Y(3)M 2 . (12) 
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Z u collects unitarity corrections corresponding to one-loop graphs with two 0(p 2 ) vertices: 
Zu = ^2 J dx dy {{d^ - g ^U} M PQ (x - y) - g iW Lp Q {x - y)}f PQ (x)f u QP {y) (13) 



P,Q ' 

-d^K PQ (x - y)f PQ {x)a QP (y) + - J r (x - y)ap Q (x)a Q p(y) 

where J, K, L, M are (renormalised) functions defined from the one-loop scalar integral with 
mesons P and Q propagating in the loop, and f 1 ^ and a = a A + a x collect source terms. 

• Za is the Wess-Zumino functional collecting anomalous contributions. 

The one-loop functional eq. (|9|) has been derived using the propagators and couplings of the 
0{p 2 ) chiral Lagrangian, and therefore it is expressed only in terms of chiral couplings: Fq and Bq, 
L{. . . [S] In particular, the Goldstone degrees of freedom have masses truncated at 0(p 2 ), denoted 

o 2 

Mp- Large fluctuations should induce significant differences between this quantity and the physical 

o 2 

mass Mp. Therefore, we want to replace Mp by Mp only when justified by physics arguments, 
since this replacement may have an important impact when comparing chiral expansions with 
experimental data. 

• The anomalous contribution Za corresponds to local couplings for vector and axial currents, 
and is not affected by our discussion. 

• For the unitarity corrections Z u , were we to consider higher and higher orders of the chiral 
expansion, we should obtain that the masses occurring in the functions J r , K, L and M r 
are physical masses, in order to get the low-mass two-particle cuts at the physical positions. 
Therefore, we write those functions with the physical masses of the Goldstone bosons. On 
the contrary, we keep the multiplying factors and a expressed in terms of parameters of 
the effective Lagrangian (m q , Bq. . .). 

• The tadpole contributions present in Z t are derived using the 0(p 2 ) contribution to the 

o 2 

Goldstone boson masses Mp. In ref. |12| . we have proposed the replacement: 

o 2 o 2 o 2 

i^log^ -> i^log^ (14) 

o 2 

We could have kept Mp everywhere in Zt, and in particular inside the logarithm. However, 
the resulting expressions are easier to deal with, and the change has only a tiny numerical 

o 2 

impact: either Mp is close to its 0(p 2 ) term and the change is trivially justified, or Mp is 
much smaller than M| and the whole tadpole contribution is very small. 

• Physical 5-matrix elements are obtained from the Green functions derived with the gener- 
ating functional by applying the LSZ reduction formula. The external legs corresponding to 
incoming and outgoing particles must be put on the mass shell. In the process, the products 
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of external momenta are translated into the well-known Mandelstam variables. These kine- 
matical relations are valid for physical masses, and we will use the latter (and not the 0(p 2 ) 

o 2 

truncated masses Mp) whenever we reexpress products of external momenta. This prescrip- 
tion is consistent with the use of physical masses in the one-loop scalar integral present in 
the unitarity term Z u . 

Following the renormalisation procedure in ref. [9], one can check easily that eq. (|14p does 
not change the renormalisation-scale dependence of LECs at one loop. Actually, the whole one- 
loop generating functional Z becomes exactly renormalisation-scale independent : when we follow 
the prescription given above, all the scale-dependent logarithms present in Zt (explictly shown in 
eq. (I10p ) and Z u (hidden in the one-loop functions M r and J r in eq. (|13|) ) are multiplied by terms 
of the same form m q Bo and thus cancel exactly. In the more usual treatment of the tadpoles [5], 
m q BQ terms are replaced by physical Goldstone masses in the one-loop generating functional (see 
sec. 8 in ref. [9]). In this case, the cancellation of the logarithms takes place only up to 0(p 4 ) and 
some higher-order logarithmic pieces of Zt have no counterpart in Z u . 

We call "bare expansion" the chiral expansion treated according to our prescription, because of 
we prefer keeping original couplings of the chiral Lagrangian to trading them for physical masses 
and decay constants. We sum up our method to obtain bare expansions of Green functions in 
Resummed xPT: 

1. Consider a subset of observables suitable for a chiral expansion, such as the linear space of 
connected QCD correlators of axial/vector currents and their derivatives away from kinematic 
singularities. 

2. Extract the bare expansion of the observables using the one-loop generating functional eq. Q: 
in Z t , replace the tadpole contributions by eq. (fT4"|) . and in Z u , use the physical masses for 
the functions J, K, L, M defined from the one-loop scalar integral. 

3. Use physical masses to reexpress scalar products of external momenta in terms of the Man- 
delstam variables. 

4. Keep track of the higher-order contributions by introducing remainders, i.e. NNLO quantities 
which have an unknown value but are assumed small enough for the chiral series to converge. 

5. Exploit algebraically the resulting relations, and never trade the couplings of the chiral La- 
grangian for observables while neglecting higher-order terms. 

The main differences from the usual treatment of three-flavour chiral series consists in the choice of 
a particular subset of observables, the distinction between physical meson masses and their 0(p 2 ) 
truncated forms, and the algebraic use of chiral expansions while keeping track of higher-order 
terms explicitly. 
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2.3 Masses and decay constants of Goldstone bosons 



The first example consists in pseudoscalar decay constants and masses. The usual expressions 
(Sec. 10 in ref. [9]) become the following bare expansions in Re^PT (similar expressions for rj can 
be found in refs. [12, 25J): 



1 K 

FlMl 
FlM\ 



F 2 Z{3) + 8(r + 2)Y(3)M 2 AL 4 + 8Y(3)M*AL 5 + F 2 e n , 
F%Z(3) + 8(r + 2)Y(3)M 2 AL 4 + 4(r + 1)Y(3)M 2 AL 5 + F 2 , 
F 2 M 2 X(3) + 16(r + 2)Y 2 (3)M*AL 6 + 16y 2 (3)M^AL 8 + F 2 M 2 d n 



r + 1 



F 2 M 2 X(3) 



(15) 
(16) 
(17) 

(18) 



+8(r + 2)(r + l)y a (3)A#AL 6 + 4(r + 1) 2 Y 2 (3)M*AL 8 + F l K M l K d K , 



We take as free parameters the Nf = 3 quark condensate and pseudoscalar decay constant expressed 
in physical units @, as well as their ratio and the quark mass ratio: 



X(3) 



2mS(3) 
~FfMf 



Z{3) 



Fl 



Y(3) 



X(3) 2mB 



Z{3) 



Ml 



m s 
m 



(19) 



We have introduced the NNLO remainders d w , e n , dx and ex, and the combinations of LECs and 
chiral logarithms: 



ALa 



AU 



ALa 



ALz 



1 



256vr 2 



log 



Ml 



128vr 2 (r-l)(r + 2) 



log 



Ml 
Ml 



+ 1 + 



2r 



1 



256vr 2 
1 



M 



Ml 



log^ + 21og^ 



Ml 



2567r 2 (r 



3 log — f- + 5 log 



M 



K 



Ml 



Ml 



1 



5127T 2 



Ml 2 Mf 
l0g ^ + 9 l0g # 



1 



512vr 2 (r + 2)(r 
1 



1) 



Mh M 2 
3 log ^+ log 



Ml 



Ml 



5127T 2 

1 



M 2 2 , M 2 
log^ + -log^ 
3 fi z 



512vr 2 (r - 1) 



3 log 



M 2 



log 



Ml 



log 



M 2 



(20) 



(21) 



(22) 



(23) 



4 In this paper, we work in the isospin symmetry limit, where m u — ma = m and the electromagnetic interaction 
is ignored. We take the following values for the masses and decay constants: Fn = 92.4 MeV, Fk/F^ — 1.22, 
= 139.6 MeV, M K = 495.7 MeV, M v = 547 MeV. 
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The values of the logarithms are only mildly dependent on r; for r = 25, 

AL 4 = L r 4 (M p ) + 0.51 • 1CT 3 , AL 5 = L r 5 {M p ) + 0.67 • 10~ 3 , (24) 

AL 6 = Lj^Mp) + 0.26 • 10~ 3 , AL 8 = Lg(M p ) + 0.20 • 10~ 3 . (25) 

Since F n , Fx, M n and Mx are accurately known, we can use these expressions to eliminate 
some of the 0(p 4 ) LECs in the chiral expansion of other observables. This is rather different from 
the usual xPT trading, since we keep explicitly higher-order terms that would have been neglected 
in the usual (perturbative) treatment of chiral series. 



From the masses and decay constants (|15l) - (jT8l) . we get the equivalent set of equations provid- 
ing some 0(p ) LECs in terms of physical masses and decay constants, r, X(3),Y(3) and NNLO 
remainders: 

^ 2 (3)AL 6 = ^-^^[l-e^-X^-d], (26) 
y 2 (3)AL 8 = l||[ e (r)+d'], (27) 



16 Ml { 

Y ^ AL * = ^T2) M^ 1 ~ V ^ ~ Z ( 3 ) ~ e ] ' ( 28 ) 
Y(3)AL 5 = I^ [r/ (r)+e']. (29) 



e(r) = 2^4, rj(r) = (% - 1 ) , (30) 



with 

1' r-l\F% 
and the following linear combinations of NNLO remainders arise: 

t + 1 / 2 \ 

d = -d n - e(r) H ) d K , d' = d - d w , (31) 

r — 1 \ r — ly 

r + 1 / 2 \ 

e = -e^ - r/(r) H ex , e' = e - e n . (32) 

r — 1 V r — 1 / 



The above identities are algebraically exact, but they are useful only as long as NNLO remain- 
ders are small. In refs. [121 125|. the size of the NNLO remainders was taken as 

d,e = 0(m 2 s ) ~ 10% , d! , e' = 0(mm s ) ~ 3% , (33) 

with the rule of thumb that NNLO corrections of size 0(m 2 s ) should not exceed (30%) 2 ~ 10% of 
the contribution to the observable while 0(m s m) terms would be less than 30% • 10% ~ 3%. We 
will propose in the next section a different but compatible way of dealing with this issue. 

In eqs. (|26|) - ()29j) . the presence of powers of Y{3), i.e., Bq, follows from the normalisation of 
the scalar and pseudoscalar sources in ref. [9]: these powers arise only for 0(p 4 ) LECs related to 
explicit chiral symmetry breaking (two powers for Lg, Lj, L%, one for L4 and L5), and are absent 
for LECs associated with purely derivative terms. 
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2.4 Prom bare expansions to Re^PT expansions 



As shown in detail in ref. [TT], plugging eqs. (|26p -(|29p i n t° the bare expansions for other observables 
corresponds to resumming the vacuum fluctuations encoded in L4 and Lq. As an illustration, we 
recall that we can exploit eqs. (f26|) - (j29|) to relate Y(3) to the chiral couplings L4 and Lq: 

F(3) = 2H-«(D-4 (M) 

[1 - V (r) -e} + - r](r) - e} 2 + k x [2AL 6 - AL 4 ] 

k = 32(r + 2)^[l-c(r)-d|. (35) 

If vacuum fluctuations are small, i.e. ALq and AL4 almost vanishing, one can treat fcx [2ALg— AL4] 
in the denominator as a small perturbation and linearise the equation as Y(3) = 1 + 0(p 2 ). This 
corresponds to the usual (iterative and perturbative) treatment of chiral series. However the factor 
k is very large (k ~ 1900 for r = 25) and values of ALq and AL4 of a few 10 -3 suffice to yield an 
important deviation of Y(3) from 1, while the linear approximation becomes inaccurate. Similar 
relations exist between X(3) and ALq, and between Z{3) and AL4 [llj . 



Using eqs. (|26j) - ([29jh we obtain the one-loop expansions of good observables in Re^PT, by us- 
ing eqs. (p?j) - P9]) and reexpressing F(3), mB , Y{3)L 4 , K(3)L 5 , Y(3) 2 L 6 , Y(3) 2 L 8 (and Y(3) 2 L 7 
through T) identities) in terms of the three parameters of interest X{3), Z(3), r and NNLO re- 
mainders. In the case of tttt and ttK scatterings, only three 0(p 4 ) LECs [L\, L2, L3) will remain. 
The square root induced by equations like eq. (|34p is a non-perturbative feature of our framework. 
It amounts to resumming (potentially) large contributions of vacuum fluctuations, encoded in the 
Zweig-rule violating LECs L4 and Lq. This feature, contrasting with the usual treatment of chiral 
series, has led us to call our framework Resummed Chiral Perturbation Theory or RexPT. 

There is a price to pay for this extension of the chiral framework in the case of large fluctuations 
of ss pairs, and the resulting competition between LO and NLO in the chiral counting: some usual 
0(p 4 ) relations cannot be exploited anymore, because of our ignorance about their convergence. 
For instance, the quark mass ratio r = m s /m (m = m u = rrid) cannot be fixed from M^/M 2 since 
we do not control the convergence of its three-flavour chiral expansion, r becomes a free parameter 
which can vary in the range: 

F K M K F 2 M 2 , , 

r '= 2 iCM7- 1 ~ 8£r£r2=2 iM- 1 ~ 3<i - (36) 

Similarly, one cannot determine LECs or combinations of LECs through ratios of observables. For 
instance, one should not use Fk/F^ to determine L5 at 0(p 4 ), because we do not know if the chiral 
expansion of Fx/F n converges at all. Finally, the agreement of the pseudoscalar spectrum with the 
Gell-Mann-Okubo formula requires a fine tuning of Lj (however, this fine tuning is also needed in 
the case of a dominant Nf = 3 quark condensate and small vacuum fluctuations |12j). 
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3 7T7T and ttK scattering amplitudes 



In this section, we are applying the Re^PT framework to two examples of Goldstone-boson scat- 
terings : 7T7T scattering, which probes the structure of QCD vacuum in the Nf = 2 chiral limit, and 
ttK scattering, which is linked with the Nf = 3 chiral limit. 

3.1 One-loop expression in Re^PT 

In the isospin symmetry limit, the low-energy irir scattering is described by a single Lorentz- 
invariant amplitude: 

A(ir a ( Pl ) + vr V) - vr c (p 3 ) + = 5 ab 5 cd A(s, t, u) + S ac 6 bd A(t, u, s) + 5 ad 5 bc A(u, t, s) (37) 

where the usual Mandelstam variables are: 

s = {p\+P2) 2 , t = (Pi-P3) 2 , u={pi-p A f, (38) 

and A is symmetric under t «-> u exchange. In a similar way, we consider the low-energy ttK 
scattering, which can be decomposed into two amplitudes according to isospin in the s-channel 
I = 3/2 and / = 1/2: 

A(7r a ( Pl ) + K\p 2 ) - 7r 6 (p 3 ) + K\ P4 )) = F* K (s,t,u) , (39) 
from which one can define two amplitudes, respectively even and odd under s <-> u exchange: 

B(s,t,u) = ^F^(s,t,u) + ^Fl /2 (s,t,u), (40) 

C(s,t,u) = -l F ^(s,t,u) + ^F^( S ,t,u). (41) 
In addition, crossing symmetry provides a relation between the two amplitudes: 

F 1 Jk(^u) = ~F*g(u,t,s) - l -F^{s,t,u) . (42) 

We can apply the prescriptions described in sec. 12.21 to determine the one-loop Re^PT expan- 
sions of A, B and C. The relevant good observables, which can be derived from Green functions 
of vector /axial currents, are F~A, F^FkF and F w FkG. 

1. We determine the one-loop bare expansions of these quantities. This can be done using 
the generating functional of Nf = 3 xPT [9], with the essential difference that we keep the 
distinction between 0(p 2 ) truncated masses and physical masses of the Goldstone bosons. 
This was performed in the case of ttK scattering in ref. [26]. A similar work can be done in 
the case of tttt scattering. The corresponding (rather lengthy) expressions are summarised in 
app. El 
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2. We use eqs. (]2T)]) - ([2"3"]) to reexpress the 0(p 4 ) LECs L4, L5, Lq, Lg in terms of r, X(3) and 
Z(3), and NNLO remainders related to n and if masses and decay constants. We denote 
with the superscript LO + NLO the resulting expressions, which include the LO and NLO 
expansions of the relevant good observables and resum the vacuum fluctuations encoded in 
L4 and Lq. 

3. To obtain the Re^PT expansions of the nir and irK scattering amplitudes, we add to the 
resulting expressions a polynomial modeling higher-order contributions : 

F ± A Re X PT = F * A LO + NLO +F 2 {sA _ M 2 )ai+F 2 {s _ SA)a2 ^ 

+(s - s A ) 2 a 3 + [(i - t A ) 2 + {u- u A ) 2 ]a± , 
F 2 F 2 B R ^ PT = F 2 F 2 K B LO + NLO + F„F K t B b l + F^F K (t-t B )b 2 (44) 

+(* - tB?h + [{s - s B ) 2 + (u- u B ) 2 ]b 4 , 
F*F$C** PT = F^ K C LO+NLO + F n F K (s-u)c 1 + (t-t B )( S -u)c 2 , (45) 

where (54, t A , ua), (s B ,t B ,u B ), (sc,tc,uc) denote the points around which we perform the 
expansion of the NNLO polynomial. The first remainder is multiplied by a constant estimating 
roughly the value of the amplitude at the expansion point (obtained from the LO chiral 
expression). The other remainders are multiplied by polynomials in the Mandelstam variables 
which vanish at the expansion point and respect the crossing properties of the amplitude. 



For our purposes, we take: 

(s A ,t A ,u A ) = (4/3M2,4/3M w 2 ,4/3M2), (46) 
(s B ,t B ,u B ) = (sc,tc,u c ) = (M 2 K + l/3Ml4/3MlM 2 K + l/3M*). (47) 



The remainders aj, 6j, Cj include only NNLO terms or higher : we expect therefore these contri- 
butions to be suppressed by 1/A^ where Ah is a typical hadronic scale [27J. On the other hand, 
the numerator may depend on the remainder considered, but the contribution to the polynomial 
must be order 0(p 6 ) in the usual chiral counting. This means that the remainders have a typical 
size of order: 

M^ 1 F 2, M 2 
ai,a<2,b 1 ,b 2 ,ci ~ — ^ , a 3 , a 4 , 63, b 4 , c 2 r • ( 48 ) 

Remainders associated with higher-order polynomials would be of order F^/A H , much suppressed 
compared to the terms considered here, and thus neglected in the following analysis. 

In the case of irir scattering, we can exploit the behaviour of the amplitude in the Nf = 2 
chiral limit in order to constrain the size of NNLO remainders further. Indeed, from Nf = 2 chiral 
perturbation theory, we know that: 

F*A(s, t, u) - F 2 (s - Ml) = 0(e 4 ) with e 2 ~ p 2 ~ m . (49) 

e counts only powers of m = m u = rrid but not those of m s . If we compare this relation with F^A 
expressed in Nf = 3 Re%PT in eqs. (|4*3"]) and ([77]) . we see that the relation ([4*9"]) implies a constrain 
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on the NNLO remainders : a\ — — (d w — e 7r )/3/(s^/M^ — 1) and a 2 — e Tl 
to m. Therefore, we can expect the remainders to exhibit the typical sizes: 



must be proportional 



d 



3(s A /M% - 1) 



,o 2 - 



M%M 2 K 
A 4 



6i,6 2 ,ci 



O3, 04,63,64,02 



A 4 h 



(50) 



According to this discussion, we take the following ranges for the direct remainders: 



FlMl 



3(s A - Ml) 



(d — e), a 2 — e G 



2MlM 2 K 
A 4 : 



6i,6 2 ,ci e 



Af/ 



A 4 h 



03,04,63,64,02 G 



2MlM 2 K 

Afj 
^ 2 Mi 



A 4 



gj4 
A 4 , 



(51) 



with Ah = 0.85 GeV. This choice for the numerical value of Ah provides a good agreement of 
our estimates with those used in refs. [U 0J E] for energy-independent quantities. In the latter 
references, NNLO remainders were taken of order 0(m 2 ) = (30%) 2 = 10% of the leading-order 
value, unless they were suppressed by one power of m and thus of order 0(mm s ) = 30% x 10% = 3%. 
According to this work, the same remainders must remain respectively of order M^/A H = 12% and 
2M 2 M^/A% = 2%. In addition, one can check that the definition and size of remainders given in 
this section can be applied to the two-point correlators related to Fp and F 2 M 2 with an expansion 
around the point of vanishing transfer momentum, leading to remainders identical to those defined 
in sec. 12.31 



3.2 Roy and Roy-Steiner equations 

The above theoretical expressions for low-energy -kit and ttK scattering must be compared to exper- 
imental information in order to extract the parameters of three-flavour chiral symmetry breaking. 
Fortunately, dispersion relations provide an appropriate framework to analyse experimental data 
and extract the low-energy behaviour of the amplitude, through Roy and Roy-Steiner equations. 

In ref. [3], Roy equations were derived and solved with experimental input on high-energy irir 
scattering. The solutions were parametrised in terms of two scattering lengths a(j and Oq. In 
refs. [3J and [I], these solutions, and some of their extensions, were exploited together with recent 
data on irir scattering in order to determine the low-energy structure of the amplitude with the 
best accuracy. Ref. [5] proposed to combine data on 5® — 5{ supplemented with a theoretical 
constraint from the scalar radius of the pion. This constraint was assessed critically in ref. [1], 
where it was proposed to avoid any reference to the scalar radius of the pion and to rely only on 
experimental data, namely data supplemented with 1 = 2 data. We follow the latter approach 
and take the results of the so-called "Global" fit, eq. (12) in ref. [1], for 7T7r-scattering data. 

In ref. [10], Roy-Steiner equations were investigated to study the ttK scattering amplitude. In 
spite of recent progress in r — » Kirv T and D — > Kirev e decays, low-energy data on ttK phase shifts 
is still lacking. However the dispersive analysis of the data in the intermediate region turned out 
to provide rather tight constrains on the low-energy irK amplitude. We use the results of ref. [TO] 
for irK scattering. 
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It is a straightforward, if tedious, exercise, to exploit the dispersive representations of the 
amplitudes A,B,C found in sec. 3 of ref. [3j and in sec. 2 of ref. [TU], and to compute them in 
subthreshold regions, where none of the dispersion integrals exhibit singularities. We checked in 
particular that our representation of the low-energy ttK amplitude was in good numerical agreement 
with the subthreshold expansion presented in sec. 6.3 in ref. |10j . 

We define the subthreshold region of interest for tttt scattering as a triangle in the Mandelstam 
plane delimited by points with (s,t,u): 

(2M 2 ,M 2 ,M 2 ), (M 2 /2,3/2M 2 ,3/2M 2 ), {M 2 J2, 3M 2 , M 2 J2) , (52) 

taking into account the symmetry of the amplitude under t — u exchange. Similarly, we define for 
ttK scattering a triangle in the Mandelstam plane with: 

(M 2 K ,2M 2 ,M 2 K ), (M 2 K ,0,M 2 K + 2M 2 ), (M 2 K + M 2 , 0, M 2 K + M 2 ) , (53) 

exploiting the symmetry or antisymmetry under s — u exchange. In each triangle, we defined 
15 points regularly spaced where we compute the scattering amplitudes. Some aspects of the 
computation, and of the correlations among the points, are covered in app. [Bj 

4 Matching in a frequentist approach 

We must match the chiral expansions of the scattering amplitudes with the experimental values 
described in the previous section. We perform this matching in a frequentist approach inspired by 
the Rfit method p3]. 

4.1 Likelihood 

We collect in a vector V our 3n observables : 

V T = [A(s u t 1 ),...A(s n ,t n ),B(s[,t[),...B(slQ,C(s'(,t , ;),...C(slO} . (54) 

Since we use the masses and decay constant identities for pions and kaons to reexpress the 
0(p 4 ) LECs in terms of Fj and F|M| through eqs. (f26l) - (f29l) . our set of theoretical parameters is: 

Parameters : r, X(3), Z(3), L[, L r 2 , L% , (55) 
Direct remainders : a%, a^, 0,3, 04, 61, &2> &4j ci> C2 , (56) 
Indirect remainders : d,d ,e,e ,dx,dz ■ (57) 

We have separated the direct remainders, attached to the bare expansions of the observables, and 
the indirect remainders, arising through the reexpression of 0(p 4 ) LECs thanks to mass and decay 
constant equalities. The latter include also the remainders dx and dz, whose expressions will be 
given in sec. 14.21 and which are required to express the paramagnetic constraints on X and Z, 
eq. ©. 
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We construct the experimental likelihood C exp , i.e. the probability of observing the data for a 
given choice of theoretical parameters T n : 

£ exp (T n ) = P(data|T n ) oc exp (-\{V th - V exp f C' 1 (V th - V exp )) /VdetC. 



(58) 

To avoid a proliferation of (purely numerical) normalisation factors of no significance for our dis- 
cussion, we use the sign oc meaning "proportional to". C is the covariance matrix between the 
experimental values V exp computed through eq. (|82j) . whereas Vth denote the theoretical values 
computed with the particular choice of T n . Since we expect strong correlations among the param- 
eters, the covariance matrix must be treated with some care, as described in app. [C] 

The theoretical likelihood £th(T n ) describes our current knowledge on the theory parameters. In 
agreement with the Rfit prescription [13], we consider that Cth(T n ) = 1 if each theoretical parameter 
lies within its allowed range described in the next section, otherwise the likelihood vanishes. 



4.2 Constraints on the theoretical parameters 

To build the theoretical likelihood, we impose a list of constraints on the theoretical parameters. 
Some constraints are fairly simple: 

• We take the following range for the ratio of quark masses r: 

FkMk 1 n f F K M K \ 2 

• Vacuum stability yields constraints on the Nf = 3 chiral order parameters: 

X(3)>0, Z(3)>0. (60) 

• We allow the three 0(p 4 ) LECs L[{M p ) , U 2 (M p ) , L 3 in the range [-F%/A%,F%/A 2 H ], i.e. lower 
than 12 • 10 -3 in absolute value. 

• The direct remainders are constrained to remain in the range given in eq. (|5ip . 

• The indirect remainders must lie in the ranges discussed in sec. 13.11 



d ,e ,d x ,e x G 



2MlM 2 K 2MlM 2 K 



d, e G 



(61) 



i.e. 3% for the first and 12% for the latter. 
A second set of constraints translates into bounds on combinations of remainders: 
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Vacuum stability for Nf = 2 chiral order parameters yields: 

X{2) > d < d max = 1 - e(r) - Y(3) 2 x L x , (62) 
Z(2)>0 e < e max = 1 - 7y(r) - Y(3) x L z , (63) 

where Lx and are small combinations of chiral logarithms denoted fx and gi in refs. [Ill 
112] . These chiral logarithms involve M% and in the Nf = 2 chiral limit, which can be 
computed through the iterative method presented in ref. [11] . 

The paramagnetic inequalities eq. © lead to: 

X(S)<X(2) - d x > d XMn = 1 - x d ™ X { ~_ d d) , (64) 

Z(3)<Z(2) <-> e z >e z , mi n^l- ~ % ■ (65) 

Z(3)(l - e) 

The ratio of order parameters Y(3) = X(3)/Z(3) = 2rnBo/M% is bound [TT] : 

1 _ fO^) _ 

y(3) < y max = 2- — -H — - (66) 
1 — ??(r) — e 



4.3 Computation of the confidence level 

Contrary to ref. [T2] which adopted a Bayesian approach to deal with tttt scattering, we follow the 
(frequentist) Rfit procedure advocated in ref. pJ3] and used for the analysis of the CKM matrix in 
ref. [28] . From the theoretical and experimental likelihoods we define the function of theoretical 
parameters 

X \T n ) = -21og£(T n ) = -21og[£ th (T n )£ exp (T n )] . (67) 

We start by computing the absolute minimum of x 2 ) letting all theoretical parameters vary 
freely : we denote X m i n a ii ^ ms vame - Then we focus on one particular theoretical parameter Tj. 
We assume that it reaches a particular value t( and compute the minimum: 

Xmin;not »(*») = min{x 2 (T n ) ; = tj} . (68) 
Then we compute the corresponding confidence level: 

V(ti) = Prob[x mi n;not ~ Xmin;alh 1] > (69) 

where Prob(c 2 , A^ f) is the routine from the CERN library providing the probability that a random 
variable having a ^-distribution with iVjof degrees of freedom assumes a value which is larger 
than c 2 . Admittedly, we are simplifying the statistical problem at hand, since we assume that the 
function x 2 (ti) nas indeed a x 2 -distribution. This should be a correct assumption if the experimental 
component is free from non Gaussian contributions and inconsistent measurements |13j . 

This method provides an upper bound on the marginal confidence level (CL) of Tj = ij for 
the optimal set of theoretical parameters : the CL value is the probability that a new series of 
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measurements will agree with the most favourable set of theoretical parameters (at T, = ij) in a 
worse way than the experimental results actually used in the analysis [J3]. The value of ti for which 
Vifi) is maximal provides an estimator of Tj : in the ideal case of very accurate data in excellent 
agreement with theoretical expectations, "P(ii) should exhibit a sharp peak indicating the "true" 
value of Tj. 

We have implemented this procedure in a program. Before turning to Goldstone boson scatter- 
ing, we checked the validity of our programs using "fake" observables. We designed observables with 
very simple chiral representations (linear or quadratic dependence on r, X(3),Z(3)) and we simu- 
lated a set of data with a certain choice of r,X(3), Z(3), adding some random noise. We plugged 
these "data" into our program and computed the confidence level for each theoretical parameter 
r, X(3),Z(3). When the chiral representation of the observables depended on this parameter, we 
obtained a function 'P(ij) showing a peak in agreement with the value used to simulate the data 
(i.e., we recovered the information contained in the data). When the chiral series for the observ- 
ables had no dependence on the parameter, the function "P(t«) was flat (i.e., we did not extract 
information absent from the data). 

5 Results 

In this section, we discuss the results obtained by matching the one-loop Re^PT expansions and 
the dispersive results on tttt and ttK scattering, relying on the frequentist approach described in 
the previous section. 

5.1 CL for order parameters and related quantities 

We have plotted the confidence level of the order parameters -^(3), Y(3) and Z(3), as well as the 
quark mass ratio r. In each case, the dashed line indicates the results obtained from tttt scattering, 
the dotted line from irK scattering, while the solid line stems from the combination of both pieces 
of information. 

If we include tttt scattering only, we see that small values of r, below 13, are disfavoured (this 
is also the case for large values of r above 25, but not at a significant level) : r > 12 at 68 % 
CL. The CL for X(3) is flat up to 0.85, where it suddenly drops, as well as that for Z(3) up to 
0.95. Y(3), which is related to -Bo an d measures the fraction of the LO contribution to M%, is 
essentially not constrained, even though values close to 2 are slightly disfavoured. If we consider 
irK scattering only, r and Y(3) are essentially not constrained. Flat CLs are observed for X(3) 
and Z(3), with a steep decrease respectively for 0.83 and 1. Finally, if we combine both pieces 
of information, intermediate values of r are clearly favoured (between 20 and 25), in agreement 
with the information contained in tttt and ttK scattering data. Low values of X(3) and Y(3) are 
preferred, whereas the CL for Z(3) peaks around 0.8. We see that the combination of the two 
data sets provides more stringent constraints on the various theoretical parameters of interest (this 
issue is discussed in more detail in app.[D]), even though these results have still a limited statistical 
significance. 
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7T7T data 


ttK data 


■kit and ttK 


Roy-Steiner 0(v A ) 


Kh, 0{p A ) 


Kh, 0(p e ) 


10 3 Li 
10 3 L 2 
10 3 L 3 


[-8.1,5.6] 
[0.2,2.4] 
nd 


[-4.4,4.1] 
nd 
nd 


[-2.1,2.2] 

[0,3.0] 
[-7.8,3.4] 


1.05 ±0.12 
1.32 ±0.03 
-4.53 ±0.14 


0.46 ±0.24 
1.49 ±0.23 
-3.18 ±0.85 


0.53 ±0.25 
0.71 ±0.27 
-2.72 ± 1.12 



Table 1: Derivative chiral couplings L T 12 z{l jL ) at [i = 0.77 GeV obtained in our approach. The 
confidence intervals correspond to a 68 % CL. "nd" means that the corresponding CL is flat over 
the whole range imposed by the theoretical likelihood, and thus the coupling is not determined. Also 
shown are results obtained assuming small vacuum fluctuations of ss pairs : ref. J775J/ analysed 
subthreshold ttK parameters from Roy-Steiner equations at order p 4 (col. 5), whereas ref. \32j 
performed fits to the KI4 form factors using chiral expansions at order p A ( col. 6) as well as p 6 
(col. 7). 



We recall that the frequentist method given here provides an upper bound on the confidence 
level (CL) for the optimal set of theoretical parameters assuming Tj = tj |13| . In the ideal case, 
we would expect the CL to peak in a very limited interval of tj, providing the "true" value of the 
corresponding theoretical parameter. In practice, we see that the chosen set of data is not accurate 
enough to provide very stringent constraints on the theoretical parameters. In such a case, the CL 
profiles can be exploited to extract a confidence interval, say at 68 % CL, i.e. a range of values so 
that the probability that the range contains the true value of the parameter is 68 %. This can be 
obtained by determining the region of parameter space where the CL curve lies above 0.32 |14j . 

From the CL profiles obtained from the combined analysis of tttt and ttK scattering, we obtain 
the following confidence intervals at 68 % CL : 

r > 14.8 , X(3) < 0.83 , Y(3) < 1.1 , 0.18 < Z(3) < 1 . [68% CL] (70) 

The values for Li,L2,L^ can also be determined in each case, and the corresponding confidence 
intervals are collected in table [H 

As a cross-check, we have also studied the case where the higher-order direct remainders are 
removed, i.e. eqs. (|43|) - (|45|) is set to zero. The corresponding CLs are sharper, but very similar in 
shape to those presented here. Therefore, the polynomial terms modeling higher order contributions 
tend to push CLs towards 1, but the qualitative features shown in figs. Q] and [2] stem mainly from 
the matching of LO and NLO terms of the Re^PT expansion to experimental information. 

The scenario mildly favoured from the matching of both irn and irK scatterings would corre- 
spond to a value of r = m s /m quite close to the canonical value r = 25. However, we emphasise 
that this agreement is rather coincidental : the latter value comes from the (perturbative) reex- 
pression of M^/M^ in terms of r, assuming that the chiral expansions of the two squared masses 
converge quickly. This assumption is not supported by our results for the quark condensate (or X), 
which exhibits some suppression when one moves from the Nf = 2 chiral limit to the Nf = 3 one, 
i.e., when m s decreases from its physical value down to zero. On the other hand, the pion decay 
constant (or Z) seems quite stable from Nf = 2 to Nf = 3, see eqs. ©-(J?]). If our results are con- 
firmed by further experimental data, we expect the usual treatment of Nf = 3 chiral expansions to 
yield unstable expansions, with significant numerical competition among terms of different orders 
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in the chiral counting. 

Such a situation is reminiscent of a scenario proposed some time ago concerning the Nf- 
dependence of the chiral structure of QCD vacuum [U [29]. The quark condensate S(iVj) and 
the decay constant F(Nf) depend on the way small eigenvalues of the Dirac operator accumulate 
around zero in the thermodynamic limit. It was conjectured that the two order parameters could 
decrease at a different rate when the number of massless flavours Nf increases : the quark con- 
densate would vanish first, followed later by the vanishing of the decay constant related to the 
restoration of chiral symmetry. The trend of our results for Nf = 3 order parameters, compared to 
Nf = 2 results, could fit such a scenario, but more data should be included in the analysis before 
we reach statistically significant CLs for the various theoretical parameters analysed here. 

5.2 Comparison with some earlier works 

5.2.1 7T7T scattering 

For 7T7T scattering, it is interesting to compare our results with ref. |12| . which shares some ideas 
and issues with the present paper. This work differs on three points from ref. [12J : we include 
ttK scattering in our analysis, we choose as observables the scattering amplitudes in subthreshold 
regions rather than the subtraction constants involved in dispersive representations, we perform 
the statistical analysis in a frequentist framework rather than a Bayesian one. 

We observe the same qualitative features in both analyses. As expected, low values of r are 
strongly disfavoured. Indeed, the analysis of currently available data on tttt scattering [1] provides 
a value of X(2), eq. Q. As illustrated in fig. 1 of ref. pQ, X(2) is related to r through the pion 
and kaon mass and decay constant identities, eqs. (|15 p - (|18p : the value of X{2) from ref. [1] favours 
the same range for the quark mass ratio as the upper plot in fig. [TJ On the other hand, we find 
that X{3) and Z{3) are only constrained through an upper bound, in numerical agreement with 
the paramagnetic inequalities X(3) < X{2) and Z(3) < Z(2). 

This agreement is particularly gratifying since the method of analysis of the present work does 
not require computing any Nf = 2 chiral order parameters or related subtraction constants like 
refs. [HE]. Moreover, one can see an improvement compared to the latter references, thanks to the 
frequentist approach chosen here. In ref. [12] . it was difficult to disentangle the effect of the data 
from that of the Bayesian priors inside a posterior p.d.f. : the so-called "reference profiles" (p.d.f.s 
from priors but no data) had to be compared to the posterior p.d.f.s (p.d.f.s from priors and data) 
to judge the impact of tttt data. In the present paper, this intricate procedure and the arbitrariness 
induced by Bayesian priors are avoided : it is clearly seen that tttt data constrains X(3) and Z(3) 
only through the values of X{2) and Z(2) and the corresponding paramagnetic upper bounds. 

5.2.2 ttK scattering 

For ttK scattering, we can compare our results with ref. [10] , where the solutions of the Roy-Steiner 
dispersion relations were used to reconstruct the amplitudes in the subthreshold region. These 
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X(3) 

Figure 1: CL profiles for r = m s /m (top) and X(3) = 2mS(3)/(F^M^) (bottom). The dashed 
line corresponds to experimental information on -kit scattering, the dotted line to ttK scattering, 
and the solid line to the combination of both sets. The two horizontal lines indicate the confidence 
intervals at 68 and 95% CL. 
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0,2 



0,4 



0,6 
Z(3) 



0,8 



Figure 2: CL profiles for Y(3) = 2mB /M% (top) and Z(3) = i^/^ (bottom). The dashed line 
corresponds to experimental information on 7T7T scattering, the dotted line to irK scattering, and 
the solid line to the combination of both sets. The two horizontal lines indicate the confidence 
intervals at 68 and 95% CL. 
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amplitudes were expanded around the point s = u, t = 0, and the coefficients of the polynomials, 
Cf- and C^-, were matched with their NLO chiral expansions in order to determine some 0(p 4 ) 
LECs. This led to a determination of L\, L2, L3 recalled in the previous section, and to a value of L4 
suggesting a significant suppression of Z{3). The value of Lq, though affected by large uncertainties, 
indicated also a suppression of X(3), stronger than that of Z{3) : 

Ref. [TO] : l I( m p) = (0-53 ± 0.39) • 10" 3 , [2L r 6 + L r 8 ](M p ) = (3.66 ± 1.52) • 1(T 3 . (71) 

Using eq. ([MD and the other results of sec. 3.1 in ref. [TO], and taking L r 8 (M p ) = (0.9±0.3)-l(T 3 [9] we 
can convert these results into the parameters of interest: following these results, X(3) would stand 
between 0.15 and 0.41, Z{3) between 0.14 and 0.92, and Y(3) between 0.44 and 1.05. Obviously, the 
low values of X(3) and Z{3) indicate that the values obtained in ref. [10], relying on the assumption 
of small vacuum fluctuations and on X{3) and Z{3) close to 1, should be reassessed relaxing this 
hypothesis. 

If our results for the combined tttt and irK data point towards a similar pattern, our analysis of 
ttK data alone provides weaker constraints than that of ref. [10] . At least two different reasons lead 
us to weaker constraints. First, we have explicitly take into account the presence of NNLO contri- 
butions which were neglected in the 0(p 4 ) analysis of ref. [TO] and which may affect significantly 
the energy-dependent part of the amplitudes. Secondly, the analysis in ref. [TO] assumes explicitly 
the smallness of vacuum fluctuations : once we drop this assumption, a smaller value of L4 (and 
thus a value of Z{3) close to 1) can be compensated by the variation of other parameters, such as 
the quark mass ratio r. These two phenomena may explain the weaker constraints observed in our 
analysis. 

5.2.3 Combined analyses 

For the combined analysis of tttt and irK scatterings, we can compare our results with refs. [201121] . 
The authors took a different approach from ours, computing NNLO chiral expansions to tttt and 
ttK scattering amplitudes, and matching with results on tttt scattering (scattering lengths) and 
irK scattering (scattering lengths and subthreshold expansion coefficients), supplemented with 
information on form factors. In agreement with the one- loop framework of ref. jS], these two- 
loop computations assume a numerical dominance of LO contributions and a quick convergence of 
Nf = 3 chiral expansions. 

In previous studies in this NNLO framework [22], the authors performed fits to pseudoscalar 
masses and decay constants [23], decays [21], and scalar form factors [TO]. In each case, the 
values of the Zweig-rule suppressed 0(p 4 ) LECs L4 and had to be fixed by hand : fits of 
similarity quality could be obtained with values of these two constants corresponding either to 
small or large vacuum fluctuations of ss pairs. For scalar form factors, values of L4 and Lq larger 
than conventionally assumed led to an improvement in the convergence of observables (fits A,B,C 
compared to fit 10, in Table 2 of ref. [19J). 

In the case of refs. [20} 121]. the authors analysed tttt and ttK scattering amplitudes in the same 
NNLO framework. The fits were not able to reproduce some observables, in particular among ttK 
subthreshold coefficients. A particular subset of subthreshold coefficients and scattering lengths led 
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to L\(Mp) ~ 0.2 • 10 -3 and Lq(M p ) negative. Such values correspond to Fq rather small compared 
to F„, with a rather unsatisfying convergence of some observables : for instance, the pion mass 
exhibits instabilities in its chiral expansion [22] . It proves difficult to draw a fully consistent picture 
for the structure of QCD vacuum in the Nf = 3 chiral limit from these results. 

Some of the problems encountered in refs. |20[ [2Tj were reassessed in ref. |30j . in particular 
the determination of NNLO LECs. Following ref. [31 j . the many 0(p 6 ) LECs are often estimated 
using resonance saturation. In ref. |30| . the specific resonance Lagrangian used in refs. [20^ 121 j was 
shown to provide values for vector-dominated LECs rather far away from the expectations based 
on irK dispersion relations, but other resonance Lagrangians failed also to reproduce these same 
results. Therefore, one may wonder whether the problems of convergence seen in |22j could stem 
from two different sources. The first one consists in the use of resonance saturation to fix 0(p & ) 
counterterms, which is already delicate in vector channels and certainly questionable in the scalar 
sector. The second one is the observed slow convergence of chiral expansions, which contradicts the 
starting assumptions of the NNLO analysis. A comparison of RexPT expansions with the NNLO 
formulae in refs. [20} 121 j should highlight how large values of the 0(p 4 ) LECs L4 and Lq might 
destabilise NNLO expansions and how the explicit resummation of vacuum fluctuations of our work 
echoes in the perturbative expansion adopted in the latter references. 

5.2.4 Lattice 

Other interesting developments are awaited from lattice simulations. The effects presented in this 
paper are related to strange sea-quarks, and can be tackled only with (2+1) dynamical fermions with 
light masses. Unfortunately, fermions with interesting chiral properties (Wilson, Ginsparg- Wilson, 
twisted-mass) [8] are still with at most two dynamical flavours. On the other hand, staggered 
fermions [33j have been exploited for simulations with (2+1) dynamical quarks, but their use is 
under much debate [Sj. The presence of the fourth root of the fermion determinant yields non- 
localities which are not understood yet : at best, recovering QCD requires taking the various 
continuum limits in a very careful way. 

A staggered version of chiral perturbation theory [35j has been developed to extract chiral LECs 
from the pseudoscalar spectrum. It attempts at reproducing the fourth-rooting of the fermion 
determinant and includes many other effects (lattice spacing, finite-volume effects, taste-breaking 
terms), leading to a number of LECs much larger than in continuum unstaggered xPT. The hope is 
that the LECs common to both theories should be identical because QCD ought to be recovered as a 
limit of lattice QCD with fourth-rooted staggered fermions. In practice [33], chiral fits to staggered 
data on the pseudoscalar spectrum must include a large number of parameters and thus are highly 
non-trivial. Mixed actions with domain-wall valence quarks and staggered sea quarks have also 
been considered to reduce the number of LECs involved in the associated chiral Lagrangian at the 
price of losing unitarity in addition to locality [36] . 

Bearing all these remarks in mind, we can focus on the following staggered values : 

Ref. [33J : r = 27.2(4) , [2L r 6 - L\]{M V ) = 0.5(1)(2) • 10~ 3 , L r A {M v ) = 0.1(2)(2) • 10~ 3 . (72) 

Combining the errors in quadrature and using eq. (I34D and the other results of sec. 3.1 in ref. [12] . 
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we can convert these results into the parameters of interest: following these lattice results, X{3) 
would stand between 0.55 and 0.95, Z{3) between 0.57 and 1.04, and Y(3) between 0.67 and 1.08, 
values which are not in striking disagreement with our results. Obviously, if the values of X(3) and 
Z(3) are on the smaller end of these ranges, i.e., if L§ and L4 are in the upper end of the range 
in ref. [33j . the assumption of small vacuum fluctuations is not correct, and the extraction of the 
LECs by the means of staggered xPT should be reassessed more carefully. 

As an alternative to such tests, which rely strongly on the usual treatment of chiral series, we 
proposed a lattice test of the size of ss vacuum fluctuations based on Re^PT in ref. [37] . We 
considered simulations with (2+1) flavours, with a strange quark mass at its physical value, but 
two u, d light quarks with identical masses rh larger than their physical values m and smaller than 
m s . The larger values of the u,d masses enhanced the impact of the vacuum fluctuations encoded 
in L4 and Lq on observables such as the masses and decay constants of pions and kaons. This led 
to a difference in the curvatures of Fp and FpMp (P = n, K) as functions of q = rh/m s , depending 
on the size of X(3) and Z{3). The effect was less pronounced in the case of Mp, obtained as the 
ratio of the two former observables, leading to a fairly linear behaviour as a function of q. 

We proposed in the same reference a test of the size of X(3) on the lattice from the pion and 
kaon spectrum, by considering the dependence on q of the ratios : 



assessed the leading finite- volume effects to conclude that large volumes (of side around 2.5 fm) 
were required to tame these effects. 

In any case, more dedicated studies on (2+1) fermions with different actions, lattice spacings 
and volumes will be required in order to draw definite conclusions from lattice simulations on the 
structure of Nf = 3 chiral vacuum. 

6 Conclusion 

Vacuum fluctuations of ss pairs can induce significant differences in the pattern of chiral symmetry 
breaking between the two conceivable chiral limits: Nf = 2 (m u = = but m s kept at its 
physical values) and Nf = 3 {m u = m& = m s = 0) . These fluctuations might lead to a paramagnetic 
suppression of the two main chiral order parameters in the Nf = 3 chiral limit, the quark condensate 
and the pseudoscalar decay constant, compared to their Nf = 2 counterparts [U HI] . Then, 
we would observe a numerical competition between leading-order (LO) and next-to-leading order 
(NLO) contributions in chiral series, through the two 0(p 4 ) LECs L4 and Lq related to the violation 
of the Zweig rule in the scalar sector. 

In order to shed light on the size of these fluctuations, we developed and modified the framework 
sketched in ref. [12] : Resummed Chiral Perturbation Theory or Re^PT. We applied it to our current 
knowledge of low-energy tttt and ttK scatterings. First, we recalled and detailed our treatment of 



Rw — 




(73) 



where F% and M% denote quantities computed on the lattice with u, d quarks of mass rh. We 
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one-loop chiral series in the case of large vacuum fluctuations: only a subset of "good" observables 
is assumed to converge globally (so that NNLO contributions are much smaller than the sum of 
LO and NLO contributions), the chiral series must be treated in a particular way to derive bare 
expansions an resum the effects of vacuum fluctuations, while NNLO remainders are introduced 
to keep track of higher-order contributions. Then, in this resummed framework, called RexPT, 
we determined the one- loop expansions for tttt and ttK scattering amplitudes. Relying on our 
current experimental knowledge, we exploited solutions of Roy and Roy-Steiner equations within 
dispersive representations to determine the values of the amplitudes in subthreshold (unphysical) 
regions where chiral expansions should converge. 

The two representations of the scattering amplitudes were matched in a frequentist approach 
(inspired by Rfit |13j). The output of this analysis are marginal CL curves, providing an upper 
bound on the confidence level (CL) for the optimal set of theoretical parameters at fixed T, = tf. the 
CL value is the probability that a new series of measurements will agree with the most favourable 
set of theoretical parameters (at Tj = tj) in a worse way than the experimental results actually 
used in the analysis |14j . 

Unfortunately, the marginal CL profiles do not provide sharp peaks and thus stringent con- 
straints on theoretical parameters at a statistically significant level. However, our results point 
towards some favoured regions of parameter space, see figs. [T] and [2j If only tttt scattering is in- 
cluded, the results obtained in earlier works [12] are recovered : small values of r are disfavoured, 
whereas -X"(3) and Z(3) are only constrained to remain below their Nf = 2 counterpart due to 
paramagnetic inequalities eq. ©. ttK scattering alone does not constrain strongly the various 
theoretical parameters, apart from setting bounds on X(3) and Z{3). The combination of the two 
pieces of information proves more interesting : the CL profile for r peaks around 23, low values of 
X(3) are preferred, whereas the CL for Z{3) exhibits a broad peak around 0.8. 

From the CL curves obtained from the combined analysis of tttt and ttK scattering, we obtain 
the following confidence intervals at 68 % CL : 

r > 14.8 , X(3) < 0.83 , Y{3) < 1.1 , 0.18 < Z{3) < 1 . [68% CL] (74) 

corresponding to the regions of parameter space where the marginal CL profiles lie above 0.32 [14] . 

The pattern of the marginal CL profiles is consistent with the scenario of significant vacuum 
fluctuations of ss pairs. It reminds one of the interesting possibility that the decrease of order 
parameters from Nf = 2 massless flavours to Nf = 3 is steeper in the case of the pseudoscalar 
decay constant F 2 (Nf) than for the quark condensate E(iVj). 

The present analysis constitutes a first attempt to analyse data with a limited statistical sig- 
nificance, and it relies strongly on the experimental results gathered on tttt and ttK scatterings. 
For tttt scattering, new results are expected from the NA48 collaboration on decays [7] and 
on the cusp in K — > 3ir [38, 39] . For ttK scattering, we hope to obtain more precise information 
from Di4 decays [3Ql IH] and r — > Kirv T decays [321 S3]. In addition, lattice studies could soon 
provide results for three light flavours with well-controlled actions in the chiral regime. These new 
high-accuracy data should shed some more light on the chiral structure of QCD vacuum, and in 
particular on its dependence on the number of massless flavours and the role played by the vacuum 
fluctuations of ss pairs. 
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A One-loop bare expansions of scattering amplitudes 



A.l 7T7T scattering amplitude 



Following the prescription in sec. 12.21 we obtain, for instance from ref. 
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where Mp denotes the leading-order pseudoscalar squared mass of the Goldstone boson P and the 
tadpole logarithm is 
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We recast the amplitude in the following form 
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where A and i3 are scale-dependent combinations of LECs : 
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which correspond to F 2 M 2 a r K7T /3 and B to F 2 f3^ n respectively, as defined in ref. [T6] . 

In the above expressions, we have replaced the bare masses by the physical masses in the 
(tadpole) logs and in the loop functions J r and M r . One can check explicitly that there is no \i 
dependence in the above expression of the amplitude : for each polynomial in s — 2M 2 ,t — 2M 2 ,u — 
2M 2 , the dependence of the LECs on the renormalisation scale fi cancels that of J r and M r . 



A. 2 tyK scattering amplitude 

We recall the expression obtained in ref. [26] for the 1 = 3/2 amplitude: 
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This expression is renormalisation-scale independent. In both tttt and ttK scatterings, the one-loop 
expressions obtained with the usual treatment of three-flavour xPT [9] are recovered if we treat 
chiral series perturbatively and neglect the (potentially large) difference between the truncated 
0{p 2 ) expressions and the physical values of the pseudoscalar masses and decay constants. 



B Computation of the amplitudes 

The amplitudes are smooth functions of the various experimental inputs. This means in particular 
that there will be significant correlations among the value of the same scattering amplitude at 
different points in the Mandelstam plane. We compute these correlations according to the following 
procedure. Let us call a k (k = 1 . . . n) the parameters describing the variations of the experimental 
inputs. To each of these parameters is attached an uncertainty (<?"&), and the correlations among 
them are encoded in a covariance matrix Dm, or equivalently, a reduced covariance matrix = 
Dki/fak&i)- We compute the mean value mj of the observables Xj's by setting all the parameters 
cifc to their central value o&: mi = Xj(ofc). Then, we vary the parameters one by one (the others 
being kept at their central value) and compute each time: 

Afc (- , a k\ a k dxi 

Af = Xi [a k + — -rrii = — x — + ... 81 
V P P da k 



28 



where p is a largish parameter (around 10), and the ellipsis denotes higher derivatives. Once this 
is done for all the parameters, we compute the covariance matrix for the observables: 

W£^« = Egif^ + - (82) 

kl kl 

The same procedure was followed in ref. [10J to determine the correlation matrix between the two 
7rif-scattering lengths. 

For the tttt scattering amplitude, we obtain the following values and errors for the amplitude at 



the limits of the subthreshold region (s,t,u) 

(2M%,M%,M%) A = 2.84 ±0.16 (83) 

(M 2 /2,3/2M 2 ,3/2M 2 ) A = -1.03 ±0.12 (84) 

(M 2 /2,3M 2 ,M 2 /2) A = -1.08 ± 0.11 (85) 

For irK scattering, we have in a similar way 

(M 2 K ,2Ml,M 2 K ) S = 4.09 ±0.64 C = (86) 

(M%,0,M% + 2M%) B = 2.96 ±0.60 C = -0.95 ± 0.03 (87) 

(M^ + M 2 ,0,M^ + M%) B = 2.61 ±0.60 C = (88) 



The zeroes of C are due to its antisymmetry under s — u exchange. Uncertainties are correlated. 



C Treatment of correlated data 



We expect strong correlations among the data points. This is reflected by the fact that the matrix 
C is nearly degenerate, and therefore cannot be inverted easily. In order to treat this problem, one 
can diagonalize H the matrix C : 



C = UDU T , D = diag(Ai,...,A n ), UU T = U T U = 1 , 
which yields the corresponding likelihood: 



C(AV) 
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/ v / det(27rC) . 



(89) 

(90) 
(91) 



Let us split the set of eigenvalues in two categories: large eigenvalues of order 1, collected in 
the diagonal matrix D, and almost vanishing eigenvalues, smaller than a cut-off and gathered in 
the diagonal matrix Dq: 

D = D + D , C = UDU T , C = UD U T . (92) 

5 In practice, we use the singular value decomposition method described in ref. [44] , which introduces two different 
rotation matrices on the left and on the right. This slight modification does not alter the procedure outlined in this 
section. 
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The eigenvalues in D$ are responsible for the near degeneracy of the matrix. In the correspond- 
ing directions, the exponential could be approximated with a Dirac distribution and would yield 
constrains on NNLO and higher-order remainders. Our approximation by a low-degree polynomial 
is expected to hold at the level of a few percent: numerically, a perfect agreement between data and 
experiment occurs already if Ay = 0(1%). Therefore, we cannot make much use of eigenvalues 
of the covariance matrix much smaller than (1%) 2 = 10~ 4 . This leads us to limit the analysis to 
the subspace where D is non-vanishing, and to define on this subspace D~ l = D^ 1 (see Ch. 2.6 in 
ref. [H] for a more detailed discussion on the relationships between singular value decomposition 
and matrix inversion). We chose to set the limit between small and large eigenvalues of order 10 -8 
(with only a very mild dependence of our results on the exact value of the cutoff). 

D Impact of combining tttt and irK data for marginal CL profiles 

As shown in sec. 15. 1} the analysis of irn data in our framework puts a lower bound on r (r > 12 
at 68 % CL), and an upper bound on X(3) and Z{3) (below 0.85 and 0.95 respectively). irK 
scattering does seem to bring only a lower bound on -X"(3). The combination of these two pieces of 
information proves much more powerful and allows one to extract CL intervals on r, X{3), Y(3) 
and Z(3). Indeed, we recall that our statistical method, inspired by the Rfit approach [13J, consists 
in the following steps : determine the absolute minimum of x 2 first, then fix a particular theoretical 
parameter a (among r, X(3), Y(3), Z{3)) and compute the corresponding relative minimum of x 2 > 
finally extract a CL (actually a P- value [H] from the difference between the two values of x 2 • This 
amounts to computing 

CL[a|data] = Max^CLfa; /x|data] (93) 

where \x collects all the remaining theoretical parameters (including Li and NNLO remainders). 
Therefore, the CL profiles obtained in our approach correspond to upper bounds on the CLs. In 
particular, it is enough that one set of theoretical parameters /j, yields the same value of x 2 as the 
absolute minimum to get CL(a) = 1. 

We have many theoretical parameters for the description of the scattering amplitudes : in 
addition to r, X(3) and Z(3), we have the 0(p 4 ) LECs £1,2,3 and (direct and indirect) NNLO 
remainders. Therefore, it is not particularly surprising that either tttt or ttK scattering alone is 
not enough to simultaneously put constrains on all these parameters : many equivalent situations 
(with identical x 2 an d thus CL) can be obtained with different sets of theoretical parameters. 
This degeneracy, in particular for the minimum x 2 1 comes from the possibility of compensating 
a variation in a by a modification of the remaining theoretical parameters /i within the allowed 
ranges, which tends to yield flat CL profiles when we consider only one amplitude. 

This underdetermination of the theoretical parameters - and the resulting degeneracy in CL 
values - is lifted once several sets of different sources are considered. In the present case, tttt data 
puts constraints on the quark mass ratio r and on some of the 0(p 4 ) derivative couplings. Because 
of these constraints, the regions of theoretical parameters with identical CLs are reduced, and thus 
the CL curves associated with ttK scattering exhibit more distinctive features. 

As an illustration of this phenomenon, we compute the CL for ttK data alone with Li,L2,L% 
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and r fixed to specific values, in order to mimic the interplay between tttt and irK data in CL 
curves. Let us set Li,L2,L^ to the values corresponding to the absolute minimum of the x 2 when 
both 7T7T and ttK data are considered 

L[(M p ) = -0.31 • 10" 3 , L r 2 (M p ) = 2.12 • 1(T 3 , L^{M p ) = —0.64 • 1CT 3 , (94) 

and let us set r to four different values r = 10, 20, 30, 35. With these theoretical parameters fixed, 
we determine CL[A(3); r, Li^^ttK] and CL[Z(3); r, Li^^lirK] profiles, which are drawn in fig. El 
We can see that the CL profiles of X(3) and Z(3) for r = 20 are very similar to the solid lines 
shown in figs. [T] and O corresponding to CL[X(3)\irir,irK] and CL[Z(3)\irir, ttK]. On the other 
hand, the curves for r = 10, 30, 35 are somewhat broader and flatter. 

The CL curves obtained for ttK scattering in figs. [Hand [2] are CL[X(3)|7riT] and CL[Z(3)|7r.?r], 
which correspond to the envelope of all the CL profiles of the form CL[X(3); r, Li^^ttK] and 
CL[Z(3); r, ^1,2,3 1 vri^] when varying r (and £1,2, 3)- This superimposition, dominated by values of 
r around 35, eventually yields the flat profiles in X(3) and Z{3) observed in figs. Q] and [2j On 
the other hand, when we combine tttt and ttK scatterings, a value of r around 20 is preferred by 
7T7T data (together with some ranges for £1,2, 3)- Therefore, the contribution from irK scattering 
to CL[X(3)|7T7r, irK] and CL[Z(3)\tt7t,itK] is close to the CL curve obtained for CL[X(3);r = 
20, _Li t 2,3 1 vrA'] and CL[Z(3);r = 20, Li^^IttK], i.e. r = 20 in fig. El This phenomenon explains 
how the information on theoretical parameters can be hidden in irK scattering, but is unveiled 
once combined with tttt scattering, lifting the degeneracy in the CLs and leading to a sharper 
determination of X(3) and Z{3). 

References 

[1] S. Descotes-Genon, L. Girlanda and J. Stern, JHEP 0001, 041 (2000) [hep-ph/9910537]. 

[2] S. Pislak et al. [BNL-E865 Collaboration], Phys. Rev. Lett. 87 (2001) 221801 
[hep-ex/0106071]; Phys. Rev. D 67 (2003) 072004 [hep-ex/0301040]. 

[3] B. Ananthanarayan, G. Colangelo, J. Gasser and H. Leutwyler, Phys. Rept. 353 (2001) 207 
[hep-ph/0005297]. 

[4] S. Descotes-Genon, N. H. Fuchs, L. Girlanda and J. Stern, Eur. Phys. J. C 24, 469 (2002) 
[hep-ph/0 112088]. 

[5] G. Colangelo, J. Gasser and H. Leutwyler, Phys. Rev. Lett. 86 (2001) 5008 [hep-ph/0 103063]. 
[6] J. Gasser and H. Leutwyler, Annals Phys. 158, 142 (1984). 

[7] L. Masetti [NA48 collaboration], Talk at the 33rd International Conference on High Energy 
Physics (ICHEP 06), hep-ex/0610071. 

[8] L. Del Debbio, L. Giusti, M. Luscher, R. Petronzio and N. Tantalo, hep-lat/0610059, 
hep-lat/0701009. 



31 





Figure 3: CL profiles for X(3) = 2mS(3) / '(F% Ml (top) and Z(3) = F$/F% (bottom), obtained 
from the combination of the experimental information on nir and irK scatterings, with Li,L2,L^ 
set to the values minimizing the complete x 2 ( see text). The three curves corresponding to four 
different (imposed) values of r = 10, 20, 30, 35. The two horizontal lines indicate the confidence 
intervals at 68 and 95% CL. 



32 



D. Becirevic, Ph. Boucaud, V. Lubicz, G. Martinelli, F. Mescia, S. Simula and C. Tarantino, 
Phys. Rev. D 74, 034501 (2006) [hep-lat/0605006]. 

Ph. Boucaud et al. [ETM Collaboration], hep-lat/0701012. 

[9] J. Gasser and H. Leutwyler, Nucl. Phys. B 250, 465 (1985). 

[10] P. Biittiker, S. Descotes-Genon and B. Moussallam, Eur. Phys. J. C 33, 409 (2004) 
[hep-ph/0310283]. 

[11] S. Descotes-Genon, L. Girlanda and J. Stern, Eur. Phys. J. C 27, 115 (2003) [hep-ph/0207337]. 

[12] S. Descotes-Genon, N. H. Fuchs, L. Girlanda and J. Stern, Eur. Phys. J. C 34, 201 (2004) 
[hep-ph/0311120]. 

[13] A. Hocker, H. Lacker, S. Laplace and F. Le Diberder, Eur. Phys. J. C 21 (2001) 225 
[hep-ph/0 104062]. 

[14] W.J. Metzger, Statistical Methods in Data Analysis, 

http://www.hef .kun.nl/ wes/stat_course/statist .pdf 

F. James, Statistical Methods in experimental physics, World Scientific 2006. 

[15] N. H. Fuchs, H. Sazdjian and J. Stern, Phys. Lett. B 238 (1990) 380, Phys. Lett. B 269 (1991) 
183, Phys. Rev. D 47 (1993) 3814 [hep-ph/9301244]. 

[16] M. Knecht, B. Moussallam, J. Stern and N. H. Fuchs, Nucl. Phys. B 457, 513 (1995) 
[hep-ph/9507319]. 

[17] B. Moussallam, Eur. Phys. J. C 14, 111 (2000) [hep-ph/9909292] ; JHEP 0008, 005 (2000) 
[hep-ph/0005245]. 

[18] S. Descotes-Genon, JHEP 0103, 002 (2001) [hep-ph/0012221]. 

[19] J. Bijnens and P. Dhonte, JHEP 0310, 061 (2003) [hep-ph/0307044] . 

[20] J. Bijnens, P. Dhonte and P. Talavera, JHEP 0401, 050 (2004) [hep-ph/0401039]. 

[21] J. Bijnens, P. Dhonte and P. Talavera, JHEP 0405, 036 (2004) [hep-ph/0404150]. 

[22] J. Bijnens, Prog. Part. Nucl. Phys. 58 (2007) 521 [hep-ph/0604043]. 

[23] G. Amoros, J. Bijnens and P. Talavera, Nucl. Phys. B 568 (2000) 319 [hep-ph/9907264]; Nucl. 
Phys. B 602 (2001) 87 [hep-ph/0 10 1127]. 

[24] J. Bijnens and P. Talavera, Nucl. Phys. B 669 (2003) 341 [hep-ph/0303103]. 

[25] S. Descotes-Genon and J. Stern, Phys. Lett. B 488, 274 (2000) [hep-ph/0007082]. 

[26] V. Bernard, N. Kaiser and U. G. Meissner, Phys. Rev. D 43 (1991) 2757. 

[27] A. Manohar and H. Georgi, Nucl. Phys. B 234 (1984) 189. 
H. Georgi, Phys. Lett. B 298 (1993) 187 [hep-ph/9207278]. 



33 



[28] J. Charles et al. [CKMfitter Group], Eur. Phys. J. C 41 (2005) 1 [hep-ph/0406184]. 

[29] J. Stern, hep-ph/9801282. 

S. Descotes-Genon and J. Stern, Phys. Rev. D 62, 054011 (2000) [hep-ph/9912234]. 

[30] K. Kampf and B. Moussallam, Eur. Phys. J. C 47 (2006) 723 [hep-ph/0604125]. 

[31] G. Ecker, J. Gasser, A. Pich and E. de Rafael, Nucl. Phys. B 321 (1989) 311. 

[32] G. Amoros, J. Bijnens and P. Talavera, Phys. Lett. B 480 (2000) 71 [hep-ph/99 12398], Nucl. 
Phys. B 585 (2000) 293 [Erratum-ibid. B 598 (2001) 665] [hep-ph/0003258]. 

[33] C. Bernard et al. [MILC Collaboration], hep-lat/0609053. 

[34] M. Creutz, hep-lat/0603020, hep-lat/0701018. 

C. Bernard, M. Golterman, Y. Shamir and S. R. Sharpe, hep-lat/0603027. 
S. Durr, PoS LAT2005 (2006) 021 [hep-lat/0509026]. 

S. R. Sharpe, PoS LAT2006 (2006) 022 [hep-lat/0610094]. 

[35] S. R. Sharpe and R. S. Van de Water, Phys. Rev. D 71 (2005) 114505 [hep-lat/0409018]. 

[36] J. W. Chen, D. O'Connell, R. S. Van de Water and A. Walker-Loud, Phys. Rev. D 73 (2006) 
074510 [hep-lat/0510024]. 

D. O'Connell, hep-lat/0609046. 

J. W. Chen, D. O'Connell and A. Walker-Loud, Phys. Rev. D 75 (2007) 054501 
[hep-lat/0611003]. 

[37] S. Descotes-Genon, Eur. Phys. J. C 40 (2005) 81 [hep-ph/0410233]. 

[38] J. R. Batley et al. [NA48/2 Collaboration], Phys. Lett. B 633 (2006) 173 [hep-ex/0511056]. 

[39] N. Cabibbo, Phys. Rev. Lett. 93 (2004) 121801 [hep-ph/0405001]. 
N. Cabibbo and G. Isidori, JHEP 0503 (2005) 021 [hep-ph/0502130]. 

G. Colangelo, J. Gasser, B. Kubis and A. Rusetsky, Phys. Lett. B 638 (2006) 187 
[hep-ph/0604084]. 

E. Gamiz, J. Prades and I. Scimemi, Eur. Phys. J. C 50 (2007) 405 [hep-ph/0602023]. 

[40] J. M. Link et al. [FOCUS Collaboration], Phys. Lett. B 607 (2005) 67 [hep-ex/0410067]. 

M. R. Shepherd et al. [CLEO Collaboration], Phys. Rev. D 74, 052001 (2006) 
[hep-ex/0606010]. 

[41] C. L. Y. Lee, M. Lu and M. B. Wise, Phys. Rev. D 46, 5040 (1992). 

B. Ananthanarayan and K. Shivaraj, Phys. Lett. B 628 (2005) 223 [hep-ph/0508116]. 

[42] R. Barate et al. [ALEPH Collaboration], Eur. Phys. J. C 11 (1999) 599 [hep-ex/9903015]. 
G. Abbiendi et al. [OPAL Collaboration], Eur. Phys. J. C 35 (2004) 437 [hep-ex/0406007]. 



34 



[43] M. Jamin, A. Pich and J. Portoles, Phys. Lett. B 640, 176 (2006) [hep-ph/0605096]. 
M. Jamin, J. A. Oiler and A. Pich, Phys. Rev. D 74, 074009 (2006) [hep-ph/0605095]. 

[44] W.H. Press, S.A. Teukolsky, W.T. Vetterling, B.P. Flannery, Numerical recipes - The art of 
scientific computing, Cambridge University Press. 



35 



